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ABSTRACT 

The results of numerical simulations of the dynamics of galactic disks, which are sub- 
merged INTO NONAXISYMMETRIC DARK MASSIVE HALO ARE DISCUSSED. GALAXY DISKS DYNAMICS 

in the nonaxisymmetric (triaxial) dark halo were investigated in detail by the high 
resolution numerical hydrodynamical methods (tvd & sph) and n-body models. the 
long-lived two arms spiral structure generates for a wide range of parameters. the 
spiral structure is global and number of turns can be 2-3 in depends of model parame- 
TERS. Morphology and kinematics of spiral structures were investigated in depends of 

THE HALO AND THE DISK PARAMETERS. THE SPIRAL STRUCTURE ROTATES SLOWLY AND THE ANGULAR 
VELOCITY VARIES IS QUASIPERIODIC. 



1. Dark massive halo and spiral structure 

The realistic models of disk galaxies require a 
massive enough halo. There are the results of nu- 
merous estimations of halo mass using various phys- 
ical approaches [SI [191 ES] • It is essential that the 
galaxy mass halo and mass of the disk are compa- 
rable and frequently, halo mass even exceeds mass 
of the disk inside of a sphere of optical radius Ropt 
{Ropt — (3 — 6) • r^i, Td — radial exponental scale of 
stellar disk [43l EOl UHl EI] ) • 

The absence of central symmetry in distribution 
of mass and potential is a feature of some galactic 
dark halo [ni[30l[3ll[32l[24l[36l[lll[6l[34]. Gener- 
ally we may consider triaxial halo, oblate or prolate 
halo [a [23]. 

Spiral structure formation in galaxies is impor- 
tant tasks of their physics. There are numerous 
mechanisms of spiral structure formation such as 
central stellar bar [35], swing amplification mecha- 
nism [4], gravitational instability [22], interactions 
with another galaxies [38] , hydrodynamical instabil- 
ities [H], induced star formation in fiocculent galax- 
ies [3]. 

In this work we discuss the possibility of spiral 
structure formation in disks which are submerged 
into nonaxisymmetric dark halo (scales a, 6, c are 



different (fig. [J) ). 

In numerical CDM models dark halo structure 
significantly differs from central symmetry [H [21 [71 
[TT| [23] . Attitude halo scales q = 6/a, s = c/a can 
reach 0.5-0.8. 

The various observational data indicate a possi- 
bility of the triaxial shape of halo. Important mark- 
ers are 

< Stellar halo objects (RR Lyrae, G,F-dwarfs, 

SDSS) [I71[30l[3ll[lQ]. 

< Tidal streams in Milky Way and in other galaxies 

(Sagittarius tidal stream, Sagittarius dSph [121 
[13], Virgo Over Density [HI [32], The Sagittar- 
ius Dwarf Galaxy [24]). 

< Polar Rings [42l [101 [36l [39l [26] . 

< Hyper Velocity Stars (HVS) [6l[l4l[l5]. 

< Galactic warps [27] . 

< Some features of galaxies structure can give the 

information about the shape DH [3ll [28l [T6l 

The most important thing for our work is nonax- 
isymmetric distribution of dark matter in the disk 



plane(fig. [T]). We will describe this feature by pa- 
rameter 

e = l--, (1) 

a 

which is limited as follows < e <0.15. We shall 
consider the effect of the shape of dark halo on the 
dynamics of the disk for £ > 0. 

2. Dynamical models of stellar and gaseous 
disks 

We have considered dynamics of stellar and 
gaseous disks in nonaxisymmetric halo separately 
for stars and gas. The basis of model of collision- 
less stellar disk is N-body dynamical system and 
Poisson equation for gravitation: 

d T' 0^ ->■ ->■ 

-T^ = -^^Fhalo{ri,t)^Fbulge{ri) , i = 1, 2, ..., A/" 

(2) 

r dr \ dr J r^ d(p'^ dz'^ 

In case of gas disk we solve the classical hydro- 
dynamics equations: 



of symmetric pseudo isothermal halo. And the sec- 
ond halo model is analog of model by Navarro, 
Frenk& White [29]: 



^ + V(.u) = 0, 


(4) 


dt 


(5) 


d fE\ Pdg 
at \g J g at 


(6) 



We used 2 numerical methods in our simulations: 

• TVD MUSCL method 

• Lagrangian method SPH with h ^ const (h is 
kernel scale) [9]. 

We also used two rigid models of the nonaxisym- 
metric halo. First is analog of potential of pseudo 
isothermal halo: 



^h{x, y, z) = AirGghoct'^ • | In(^) 



arctg(0 1 1 + ^^ 



(7) 



where i — \\ ~2 ^ — 2 ^ — 2' ^"^ ^^^ ^^^^ ^^ equal 

y ^x ^y ^z 

values of the scales a = 6 = c, we have a usual model 



gh{x,y,z) = 



QhO 



c(i+e)' 



(8) 



That kinds of potentials reproduce flat rotation 
curve with many values of their parameters. 

We took into account solid-body rotation of these 
potentials with angular velocity ft halo- Halo ro- 
tation is slower than rotation of disk on periph- 
ery (Qdisk)' 

^halo < ^disk{Ropt) • (9) 

3. Stellar disk dynamics 

Our simulations start with axisymmetrical halo 
e{t = 0) = 0. Then halo became nonaxisymmetric 
during some time interval r, which are 1-3 periods 
of disk rotation (fig. [2j). 

In the fig. [6] radial profiles of surface density a, 
rotation curve F, dispersion of radial velocities of 
particles c^, disk vertical scale are shown for the 
stellar disk. These distributions ensure the gravita- 
tional stability of the disk. Our initial radial distri- 
bution is stable because we use Toomre parameter 
Qt = Cr/{3.36Ga/de) exceeding 2. All parameters 
do not change with time in case of axisymmetrical 
halo (fig. [ID. 

The situation changes in the case of nonaxisym- 
metric halo. Fig. [3] show the distributions of the 
logarithm of the surface density in the disk at dif- 
ferent times. 

The amplitudes of Fourier-harmonics tti = 2 at 
various radii are more complicated. We observe 
quasiperiodic changes of Fourue-amplitudes for two- 
arm wave (fig. [4]). As the result a rotation velocity 
of spiral pattern is non-stationary (fig. [5j). 

Angular velocity of rotation Qp is determined by 
spiral wave phase O: 



" ~ dt • 



(10) 



Two stages of disk dynamics recognized well. 
Slow rotation Qpi related with dynamic of cen- 
tral region and fast rotation Qp2 corresponded with 
outer disk part. 

4. Gaseous disk dynamics 

In fig. [8] the typical evolution of a gas disk is 
shown. It is the logarithm of the surface density 



in different times. The formation of two- arm trail- 
ing spiral structure is observed. For typical condi- 
tions perturbations grow to nonlinear stage during 
1-3 circulation and form system of shocks. Analy- 
sis of gas flow structure demonstrates it's difficult 
character. We also observe different shock waves 
and regions with share flow. Morphology of spiral 
structures depends on rotation curve F(r), distri- 
bution of density in dark halo ^^(r, (/?), radial dis- 
tribution of stellar disk density (j{r) and the sound 
speed. Specifically in central part of the disk often 
liding spirals forms. On periphery this spirals al- 
ways become lagging spirals. Also liding spirals can 
produce 0-like structures (some times nested 0-like 
structures, for example see fig (10] with M = 10). 

Radial profiles of surface density are drawn on 
fig. [9] for a fixed azimuthal angle. The gas disk re- 
sponds even to 1 percent of halo non-axisymmetry. 
We see a small amplitude here. But for 5 percents 
of non-axisymmetry and more shock waves appear. 

Effects of gas temperature on the gas spiral pat- 
tern shown in the fig. [TOl Different values of the 
Mach number are used. The circular rotation curves 
is fixed in these models. The pitch angle of spirals 
decreases in more cold disk. 

Kinematics of spiral pattern in the gas disk is sim- 
ilar to the kinematics of the stellar disk submerged 
into nonaxisymmetric halo. We observe character- 
istic quasiperiodic curves like in case of stellar disk. 

Effect of selfgravitation increases the efficiency of 
spiral structure formation by the nonaxisymmetric 
halo. We simulate only gravitational stable disks. 
Morphology of waves changing and its amplitude 
grows. 

We also demonstrate formation of similar spi- 
ral structures in 2d and 3d numerical simulations 
(fig ITT)) . Vertical distribution of volume density 
demonstrate vertical structure of density waves 
along radius (figlT2]). Density waves marked by red 
arrows and dots at the disk plane ( fig JTTI down) and 
also in vertical distribution at the same moment. 



4. The spiral pattern rotates slowly. Position of 
corotation radius is at the outer part of the 
disk. 

5. A possibility to apply the mechanism consider 
for the formation of inner spiral structure re- 
mains open. It is also possible that the non- 
axisymmetric halo may explain the outer spi- 
ral patterns often observed in HI line and/or 
in UV (GALEX). 

We would like to thank A.Zasov and V.Korchagin 
for their useful discussions. Numerical computa- 
tions were carried out on CKIF-MGU "Chebyshov" 
at the Moscow State University, Moscow, Russia. 
This work has been supported by RFBR A/'09-02- 
97021. 



5. Discussion and Conclusions 

1. Possibility of 2-arms spiral structure formation 
in the stellar and gas components was shown 
in case of nonaxisymmetric dark massive halo. 
Properties of spiral pattern depends on param- 
eters of disk and dark halo. 



2. Spiral patterns changes in time. 

3. In general the spiral patterns have a similar 
morphology for 2D and 3D models. 
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Figure 1: Galaxy disk submerged into triaxial halo 
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Figure 2: Nonaxissymmetric index as function of time 
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Figure 3: The logarithm of surface density in stehar disk model 




Figure 4: Evolution of Fourier harmonics amplitude for two-arms mode in different radii Value r = 1 
corresponds Avd. 
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Figure 5: Evolution of Fourier harmonics phase O for two-arms mode on fixed radius 
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Figure 6: Initial distributions of the basic parameters of the stellar disk 
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Dispersion of radial velocity 
on different radii 
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Figure 7: Evolution of stellar disk parameters in model with the axisymmetrical halo are shown. There 
is no growth of perturbations. 
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Figure 8: Evolution of gaseous disk submerged into nonaxissymmetric halo 
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Figure 9: Radial distributions of surface density along azimuth angle (/? = 0° for models with different e. 
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Figure 10: Distributions of surface density in gaseous disks for different values of Mach number M 
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Figure 11: It is a comparison of spiral patterns in 2D and 3D models. 
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Figure 12: The vertical structure of gaseous disk submerged into triaxial dark halo. 
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